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Larger and faster: revised properties and a shorter orbital period 
for the WASP-57 planetary system from a pro-am collaboration"^ 
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ABSTRACT 

Transits in the WASP-57 planetary system have been found to occur half an hour earlier than 
expected. We present ten transit light curves from amateur telescopes, on which this discovery 
was based, thirteen transit light curves from professional facilities which confirm and refine 
this finding, and high-resolution imaging which show no evidence for nearby companions. 
We use these data to determine a new and precise orbital ephemeris, and measure the phys¬ 
ical properties of the system. Our revised orbital period is 4.5 s shorter than found from the 
discovery data alone, which explains the early occurrence of the transits. We also find both 
the star and planet to be larger and less massive than previously thought. The measured mass 
and radius of the planet are now consistent with theoretical models of gas giants containing 
no heavy-element core, as expected for the sub-solar metallicity of the host star. Two transits 
were observed simultaneously in four passbands. We use the resulting light curves to measure 
the planet’s radius as a function of wavelength, finding that our data are sufficient in principle 
but not in practise to constrain its atmospheric properties. We conclude with a discussion of 
© 0000 RAS the current and future status of transmission photometry studies for probing the atmospheres 

of gas-giant transiting planets. 
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1 INTRODUCTION 


Although the first tra nsiting extrasolar planet (TEP) was only dis - 
covered in late 1999 dHeruv et al.ll200d;I Charbon neau et alj|2000h . 
and the second as recently as 2003 ( Konacki et alJ2003l) . the num¬ 
ber currently known has already exceeded 120(o The great ma¬ 
jority of those are small objects observed using the NASA Kepler 
satellite; validation of the planetary nature of these bodies has been 
grea tly helped hy the ir occurrence in systems of multiple planets 
(see lRowe et alJl^l4l hut detailed studies are difficult due to their 
small size and long orbital periods (Porbs). 

A significant number (231 as of 2015/07/21) of the known 
TEPs are hot Jupiters, adopting a definition of mass Mb > 
O.SMjup and Porb < 10 d. These are much better suited to char¬ 
acterisation with existing facilities, as their relatively large masses 
and radii, short orbital periods, and bright host stars make photo¬ 
metric and spectroscopic observations easier and more productive. 
Perhaps the single most important observable property of a planet 
is its orbital period; the period distributions of exoplanets provide 
an insight i nto the mechanisms gov e rning their formation and evo¬ 
lution (e.g. Mordasini et al.ll20093 l5 [Benftez-Llambav et al.ll201ll ; 
iHellier et alJl2012h . and a precise value is mandatory for perform¬ 
ing follow-up observations. 

In this work we study the WASP-57 system, who se planetary 
natur e was discovered by the SuperWASP consortium dFaedi et al.l 
I 201 I hereafter F13). WASP-57 contains a star slightly cooler 
and less massive than the Sun (Tefi = 5600 ± 100 K, Ma = 
0.89 ± 0.07 Mq) orbited by a planet which is the same size but 
less massive than Jupiter (Mb = 0.64 ± 0.06 Mjup, Rh = 
1.05 ± 0.05 Rjup). The moderately different properties found by 
F13 placed the planet at the lower edge of the distribution of gas 
giant TEPs in the mass-radius diagram, making it a good candi¬ 
date for hosting a heavy-element core despite the subsolar metal- 
licity of the host star ([Fe/H] = —0.25 ± 0.10). The analysis b y 
F13 was based on SuperWASP photometry jPollacco et al.ll200^ . 
radial velocities from CORALIE spectra, plus two complete and 
one partial transit light curves from the Euler and TRAPPIST tele¬ 
scopes at ESO La Silla, Chile. No further work on this system has 
been published. 

Early in the 2014 observing season a group of amateur as¬ 
tronomers noticed that the transits of WASP-57 were occurring half 
an hour earlier than expected, which is a significant fraction of the 
2.3 hr total transit duration. This was immediately confirmed using 
a transit of WASP-57 which had been serendipitously observed on 
2014/05/18 using the Danish 1.5 m telescope, at La Silla. A second 
transit observation was scheduled on 2014/06/24 with the Danish 
telescope, the ESO 2.2 m telescope and GROND imager, and the 
immediately following transit with the amateur observers in Eu¬ 
rope, allowing its early arrival to be reconfirmed. We also possess 
high-precision light curves of WASP-57 obtained in the 2012 sea¬ 
son, before the imprecision of the original orbital ephemeris be¬ 
came apparent. In this work we present all the data we have ob¬ 
tained for WASP-57, produce a revised ephemeris which can be 


* Based on data collected by MiNDSTEp with the Danish 1.54 m tele¬ 
scope, and data collected with GROND on the MPG 2.2 m telescope, both 
located at ESO La Silla. 

^ See TEPCat (Transiting Extrasolar Planet Catalogue; ISouthwo^l201 ih 
at: http : / /www .astro . keele .ac.uk/jkt/tepcat/ 


used for follow-up observations in future, measure the physical 
properties of the system to high precision, and search for variations 
of the measured planetary radius with wavelength 


2 OBSERVATIONS AND DATA REDUCTION 

A total of ten transit light curves were obtained by LB, GC, CL, 
AM and GM using telescopes of apertures between 180 mm and 
300 mm, sited in Italy. Further details of the observational setup 
and numbers of datapoints are given in Table[T] 

Two complete transits of WASP-57 were observed using the 
1.54 m Dan ish Telescope a nd DFOSC instrument at ESO La Silla, 
Chile (see IPominik et aP 1201 Oti . on the dates 2014/05/18 and 
2014/06/24. DFOSC has a plate scale of 0.39" pixel” ^ and a 2048^ 
pixel CCD, giving a field of view of 13.7' x 13.7'. We windowed the 
CCD down to 1100 x 900 and 1045 x 920 pixels to shorten the dead 
time between exposures, resulting in images containing WASP- 
57 and six decent comparison stars. Both transits were obtained 
through a Bessell R filter. The instrument was defocussed in order 
to improve the efficie ncy of the observations, and to combat time- 
correlated noise (see ISouthworth et al.] 120090 . The telescope was 
autoguided to limit pointing drifts to less than five pixels over each 
observing sequence. An observing log is given in Table|^and the 
light curves are plotted individually in Fig.[T] 

The transit on 2014/06/24 observed with DFOSC was also 
monitored using GROND dOreiner et alj bOOsh mounted on the 
MPG 2.2 m telescope at La Silla, Chile. GROND was used to obtain 
light curves simultaneously in four passbands, which approximate 
the SDSS g, r, i and z bands. The small field of view of this in¬ 
strument (5.4'x5.4' at a plate scale of 0.158" pixel”^) meant that 
fewer comparison stars were available. The instrument was defo¬ 
cussed and the telescope was autoguided throughout the observing 
sequence. Poor weather conditions (high wind) forced closure of 
the telescope before the transit finished, so the light curves have 
only partial coverage of the transit (see Tablej^and Fig.[TJ. 

We now turn to observations obtained prior to the concep¬ 
tion of the current work. We observed WASP-57 on the night of 
2012/05/10 using the BUSCA instrument on the 2.2 m telescope 
at Calar Alto Astronomical Observatory. BUSCA is capable of 
observing simultaneously in four passbands, for which we chose 
Gunn u, g, r and 2 . The motivation for these choices, and a de¬ 
tailed discussion on the use of BUSCA for planet ary transit obser¬ 
vations, can be found in ISouthworth et al] d2012h . All four CCDs 
on BUSCA have a plate scale of 0.176" pixel” but were operated 
with 2x2 binning. Whilst the full field of view of 12x 12 arcmin 
was accessible in the u band, the available field in the g, r and 
2 bands was vignetted into a circle of diameter of approximately 
6'. The instrument was defocussed and the telescope was autogu¬ 
ided throughout the observations (TableO. The u-band data are of 
insufficient quality for full modelling but can be used to obtain a 
time of minimum and to check for a possible variation of measured 
planetary radius with wavelength. 

One transit of WASP-57 was observed on 2012/04/01 with 
EulerCam, using the same methods as for the EulerCam transit 
in F13. EulerCam is a CCD imager mounted on the L2m Euler- 
Swiss telescope. La Silla, with a field of view of 14.7'x 14.7' 
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Table 1. Instrumental setup for the amateur observations. A^obs is the number of observations. 


Date 

Observer 

Telescope 

CCD 

Eilter 

Nobs 

2014/05/24 

C. Lopresti 

180 mm Maksutov-Newton 

SBIG STIOXME 

red 

39 

2014/06/10 

G. Corfini 

200 mm aperture, 800 mm focal length 

SBIG STT-1603 

deal* 

87 

2014/06/10 

C. Lopresti 

300 mm aperture, 1500 mm focal length 

SBIG STIOXME 

red 

42 

2014/06/10 

A. Marchini 

300 mm Zen Maksutov-Cassegrain 

STL-6303 

Cousins R 

58 

2014/06/27 

L. Barbieri 

300 mm aperture, 3000 mm focal length 

SBIG ST9 

clear 

40 

2014/06/27 

G. Corfini 

200 mm aperture, 800 mm focal length 

SBIG STT-1603 

deal* 

34 

2014/06/27 

C. Lopresti 

180 mm Maksutov-Newton 

SBIG STIOXME 

red 

16 

2014/06/27 

C. Lopresti 

300 mm aperture, 1500 mm focal length 

SBIG STIOXME 

red 

58 

2014/06/27 

G. Marino 

250 mm aperture, 1200 mm focal length 

SBIG ST7-XME 

clear 

64 

2015/05/28 

C. Lopresti 

180 mm Maksutov-Newton 

SBIG STIOXME 

red 

30 


Table 2. Log of the observations obtained from professional telescopes. A^obs is the number of observations, Texp is the exposure time, T^ead is the dead 
time between exposures, ‘Moon ilium.’ is the fractional illumination of the Moon at the midpoint of the transit, and A^poly is the order of the polynomial fitted 
to the out-of-transit data. The aperture radii are target aperture, inner sky and outer sky, respectively. The ‘bb’ filter is a blue-blocking filter. 


Instrument 

Date of 

first obs 

Start time 
(UT) 

End time 
(UT) 

Nobs 

Texp 

(s) 

^dead 

(s) 

Filter 

Airmass 

Moon 

ilium. 

Aperture radii 
(pixels) 

Npoly 

Scatter 

(mmag) 

TRAPPIST 

2012/03/15 

04:34 

08:20 

367 

20 

10 

bb 

1.831.121.13 

0.463 

11.3 22.7 36.3 

1 

2.25 

TRAPPIST 

2012/04/01 

04:40 

09:36 

701 

15 

7 

bb 

1.361.12-:. 1.45 

0.649 

13.5 19.3 30.9 

1 

3.49 

Euler 

2012/04/01 

05:20 

09:35 

212 

50-180 

16 

r 

1.231.121.44 

0.649 

24 



0 

0.98 

BUSCA 

2012/05/10 

23:16 

02:49 

82 

120 

35 

u 

1.31 1.29-:. 1.86 

0.665 

15 

25 

45 

1 

3.46 

BUSCA 

2012/05/10 

23:16 

03:37 

100 

120 

35 

9 

1.311.29-:. 2.46 

0.665 

17 

27 

50 

1 

1.26 

BUSCA 

2012/05/10 

23:16 

03:34 

99 

120 

35 

r 

1.31 1.29-:. 2.44 

0.665 

18 

28 

50 

1 

0.82 

BUSCA 

2012/05/10 

23:16 

03:37 

98 

120 

35 

z 

1.31-> 1.29->2.46 

0.665 

18 

28 

50 

1 

1.58 

DFOSC 

2014/05/18 

01:58 

06:51 

162 

100 

9 

R 

1.32-> 1.12-> 1.70 

0.761 

12 

20 

40 

1 

0.74 

DFOSC 

2014/06/24 

23:33 

04:15 

155 

100 

8 

R 

1.27-> 1.12-> 1.53 

0.050 

16 

26 

45 

1 

0.71 

GROND 

2014/06/24 

23:38 

02:34 

53 

90-110 

54 

9 

1.25-> 1.12-> 1.17 

0.050 

40 

60 

90 

1 

0.89 

GROND 

2014/06/24 

23:38 

02:34 

53 

90-110 

54 

r 

1.25-> 1.12-> 1.17 

0.050 

40 

60 

90 

1 

0.58 

GROND 

2014/06/24 

23:38 

02:34 

53 

90-110 

54 

i 

1.25-> 1.12-> 1.17 

0.050 

32 

55 

75 

1 

0.58 

GROND 

2014/06/24 

23:38 

02:34 

53 

90-110 

54 

z 

1.25-> 1.12-> 1.17 

0.050 

30 

50 

80 

1 

1.02 


at 0.23^^pixel“^. We obtained 212 images through a Gunn r fil¬ 
ter, without applying a defocus to the instrument. Further details 
on EulerCam and the data reduction proceduce can be found in 
iLendl et alJbOlA . 

Two transits of WASP-57 were observed on 2012/03/15 and 
2012/04/01 using the 0.6 m TRAPPIST roboti c telescope located 
at La Silla jGillon et alj|201 ihijehin embOl ih . The 2kx2k CCD 
was thermoelectrically cooled and yielded a field of view of 
22'x22' at 0.65" pixel“^. Images were obtained with a slight de¬ 
focus and through a blue-hlocking filtel^ that has a transmittance 
greater than 90% from 500 nm to heyond 1000 nm. 


2.1 Data reduction 

The data from the amateur telescopes were all reduced using 
Maxim diB In each case the science images were calibrated using 
dark and flat-field frames. 

The data from DFOSC, CROND and BUSCA were re¬ 
duced usin g aperture photometry as im plemented in the DE¬ 
POT code dSouthworth et alj 1 20091 l2014ll . which relies on the 


^ http://WWW.astrodon.com/products/fiIters/exoplanet/ 
® http: / /WWW. cyanogen. com/maximjnain. php 


IDiB/ASTROLIe|3 implementation of DAOPHOT dStetsonl fl987h . 
Master bias and flat-field images were constructed but generally 
found to have an insignificant effect on the quality of the photom¬ 
etry. Image motion was tracked by cross-correlating individual im¬ 
ages with a reference image. 

We obtained photometry on the instrumental system using 
software apertures of a range of sizes, and retained those which 
gave light curves with the smallest scatter (Tahlej^. We found that 
the choice of aperture size affects the scatter hut not the shape of 
the transit in the final light curve. The instrumental magnitudes 
were then transformed to differential-magnitude light curves, nor¬ 
malised to zero magnitude outside transit using first-order poly¬ 
nomials (Table|3 fitted to the out-of-transit data. The differential 
magnitudes are relative to a weighted ensemble of typically five 
(DFOSC) or two to four (CROND) comparison stars. The compar¬ 
ison star weights and polynomial coefficients were simultaneously 
optimised to minimise the scatter in the out-of-transit data. 

Finally, the timestamps for the datapoints were converted 
to the BJD(TDB) timescale l lEastman et al.ll20l^ . We performed 
manual time checks for several images obtained with DFOSC 


^ The acronym IDL stands for Interactive Data Language and is a 
trademark of ITT Visual Information Solutions. For further details see: 
http://www.exelisvis.com/ProductsServices/IDL.aspx 
® The ASTROLIB subroutine library is distributed by NASA. For further 
details see: http: //idlastro . gsf c . nasa . gov/. 
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Figure 1. All light curves from professional facilities presented in this work, grouped and colour-coded according to the telescope used. The inshaiment and 
hlter are labelled individually for each light curve. The second light curve from TRAPPIST has a discontinuity shortly before the midpoint of the transit due 
to a meridian-flip. This is indicated using a vertical black line. 


and verified that the FITS file timestamps are on the UTC sys¬ 
tem to within a few seconds. In recent work on the WASP-103 
system we found that the timestamps from DFOSC and GROND 
agree to within a few s econds, supporting the reliability of both 
dSouth worth et al.l2015h . The reduced data are given in Table[3and 
will be lodged with the CD^l 

The data from EulerCam were re duced using apert ure pho¬ 
tometry following the methods given bv iLendl et al.l ( l2()l^ . Differ¬ 
ential aperture photometry was also used on the TRAPPIST data, 
using carefully selected extraction apertures and reference stars. 
For more details on the TRAPPIST data reduction procedures, see 
e.g. lGillon et al.l droTsh . The transit on 2012/04/01 was obtained in 
two sequences separated by a meridian flip, and the two sets of data 
were reduced independently. 

2.2 High-resolution imaging 

We obtained several images of WASP-57 with DFOSC in sharp 
focus, allowing us to check for the presence of faint nearby 
stars whose light migh t act to decrease the observed transit depth 
foaemgen et al.ll2009h . The closest stars we found on any image 

® http://vizier.u-strasbg.fr/ 


are much fainter than WASP-57, and are over 45” distant, so are 
too far away to affect our photometry. 

We also obtained a high-resolution image of WASP-57 us- 
ing the Lucky Imager (LI) mounted on the Danish telescope (see 
ISkottfelt et al.ll20iAl2015h . The LI uses an Andor 512x512 pixel 
electron-multiplying CCD, with a pixel scale of 0.09” pixeD*^ giv¬ 
ing a field of view of 45” x 45”. The data were reduced using a ded¬ 
icated pipeline and the 2% of images with the smallest point spread 
function (PSF) were shifted and added to yield a combined image 
whose PSF is smaller than the seeing limit. A long-pass dichroic 
was used, resulting in a response function which approximates that 
of SDSS i+z. An overall exposure time of 600 s corresponds to 
an effective exposure time of 12 s for the best 2% of the images 
(Fig.lHl. The FWHM of the PSF is 4.0x4.2 pixels (0.36” x 0.38”). 

Two faint stars were detected on the LI image, at angular dis¬ 
tances of 10.99 ± 0.05” and 21.56 ± 0.07” from WASP-57, and 
fainter by 8.7 ± 0.7 mag and 8.2 ± 0.4 mag. Neither of these stars 
is sufficiently bright and close to WASP-57 to affect the analysis 
presented in the current work. We assessed the limiting magnitude 
of the LI image by placing a box with sides equal to the FWHM of 
the star on each pixel on the image. The standard deviation of the 
counts within each box was calculated, and a 3 ct detection threshold 
was generated. The mean detection threshold at a given radius from 
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Table 3. Sample of the data presented in this work (the first and last dat- 
apoints of each light curve). The full dataset will be made available at the 
CDS. 


Instrument 

Filter 

BJD(TDB) 

Diff. mag. 

Uncertainty 

TRAPPIST 

bb 

2456001.690740 

0.0014100 

0.0033979 

TRAPPIST 

bb 

2456001.847510 

0.0018100 

0.0019776 

TRAPPIST 

bb 

2456018.694550 

0.0062000 

0.0041227 

TRAPPIST 

bb 

2456018.900580 

-0.0042000 

0.0041094 

Euler 

r 

2456018.722460 

0.0015300 

0.0014096 

Euler 

r 

2456018.899430 

0.0006000 

0.0010572 

BUSCA 

u 

2456058.476778 

0.0011721 

0.0027719 

BUSCA 

u 

2456058.624676 

-0.0073014 

0.0050536 

BUSCA 

9 

2456058.476778 

0.0008973 

0.0010736 

BUSCA 

9 

2456058.657776 

-0.0006933 

0.0015483 

BUSCA 

r 

2456058.476778 

0.0004552 

0.0007219 

BUSCA 

r 

2456058.655976 

0.0005960 

0.0009239 

BUSCA 

z 

2456058.476778 

0.0018345 

0.0015082 

BUSCA 

z 

2456058.657776 

0.0003519 

0.0017916 

DFOSC 

R 

2456796.588810 

0.0004779 

0.0006752 

DFOSC 

R 

2456796.792254 

0.0006067 

0.0008685 

DFOSC 

R 

2456833.486396 

-0.0006900 

0.0007120 

DFOSC 

R 

2456833.683083 

-0.0007523 

0.0007212 

GROND 

9 

2456833.490392 

0.0000938 

0.0008559 

GROND 

9 

2456833.611846 

0.0097990 

0.0009372 

GROND 

r 

2456833.490392 

0.0000128 

0.0005593 

GROND 

r 

2456833.611846 

0.0114375 

0.0006182 

GROND 

i 

2456833.490392 

0.0002314 

0.0005610 

GROND 

i 

2456833.611846 

0.0107411 

0.0009272 

GROND 

z 

2456833.490392 

0.0010646 

0.0010125 

GROND 

z 

2456833.611846 

0.0082877 

0.0010501 


the target star was then converted to a relative magnitude. Further 
details of the detection and reduction methods are given in Evans 
et al. (2015, in preparation). The contrast curve is shown in Fig.[3 


3 TRANSIT TIMING ANALYSIS 

The issue which brought WASP-57 to our attention was the offset 
between the predicted and actual times of transit. We have there¬ 
fore obtained as many measured times of mid-transit as possible. 
We first modelled the two DFOSC transits individually using the 
JKTEBOP code (see below), as these are the two light curves which 
have full coverage of a transit with a low scatter in the data. We 
scaled the errorbars for each light curve to yield a reduced of 
xJ = 1-0 versus the fitted model. This step is necessary because 
the uncertainties from the APER algorithm tend to be moderately 
too small. 

We then modelled the light curves from the amateur telescopes 
with JKTEBOP but fitting for only the time of mid-transit and the 
out-of-transit brightness of the system. The other photometric pa¬ 
rameters were fixed to the best-fitting values from the two DFOSC 
light curves. The uncertainties in the transit times were multiplied 
by to account for the underestimated observational errors 

in most of the datasets. We performed the same process on the 
GROND data, except this time we fitted the out-of-transit bright¬ 
ness as a linear function of time rather than just a constant offset 
from zero differential magnitude. 

We then turned to published data. The discovery paper of 
WASP-57 (FI3) contains two light curves observed with TRAP- 
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Figure 2. High-resolution Lucky Image of the field around WASP-57. The 
upper panel has a linear flux scale for context and the lower panel has a 
logarithmic flux scale to enhance the visibility of any faint stars. Each image 
covers 8” X 8” centred on WASP-57. A bar of length 1” is superimposed 
in the bottom-right of each image. The image is a sum of the best 2% of the 
original images. 


PIST and one with the Euler telescope. We fitted these as above, 
with the photometric parameters fitted for the two light curves with 
complete transit coverage and fixed for the TRAPPIST light curve 
which only contains the second half of a transit. We also included 
one transit time ob tained by U. Dittler a nd lodged on the Exoplanet 
Transit Databas^ jPoddanv et al.ll20l3) . 

The SuperWASP data which triggered the discovery of the 
planetary nature comprise approximately 30000 datapoints ob¬ 
tained during the 2008, 2009 and 2010 observing seasons. These 
data were obtained and separated into individual seasons, then fit- 


^ The Exoplanet Transit Database (ETD) can be found at: 
http: //var2 .astro . cz/ETD/credit. php, see also TRESCA at: 
http://var2.astro.cz/EN/tresca/index.php. 
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Figure 4. Plot of the residuals of the timings of mid-transit for WASP-57 versus a linear ephemeris (see Table|4). The points are colour-coded according to 
their source: black for the WASP data, green for the amateur timings in the current work, blue for DFOSC, red for GROND, off-yellow for BUSCA, and grey 
for the TRAPPIST and Euler telescopes. The dotted lines show the la uncertainty in the ephemeris as a function of cycle number. 



Figure 3. Constrast curve giving the limiting magnitude of the LI obser¬ 
vation as a function of angular distance from WASP-57 (dark red circles 
connected by a red line). The closest detected star is shown as a blue data- 
point. 

ted with JKTEBOP in the same way as for the data obtained using 
amateur telescopes. The resulting season-averaged times of mini¬ 
mum are consistent with the linear ephemeris found below, but are 
of low precision. WASP-57 A, at F = 13.04, is comparatively faint 
for the SuperWASP telescopes so suffers from a large scatter in its 
light curve. We included these times of minimum light in the fol¬ 
lowing analysis, but note that they do not have a significant effect 
on the results. 

All times of mid-transit were then fitted with a straight line 
versus cycle number to determine a new linear orbital ephemeris. 
Table|4] gives all transit times plus their residual versus the fitted 
ephemeris. We chose the reference epoch to be that for our BUSCA 
observations, in order to limit the covariance between the reference 
time of minimum and the orbital period. The resulting ephemeris is 

To = BJD(TDB) 2 456 058.54910(16) + 2.83891856(81) x E 

where E gives the cycle count versus the reference epoch and the 
bracketed quantities indicate the uncertainty in the final digit of the 
preceding number. This orbital period is 4.5 s (24(j) smaller than 
the value of 2.838971 (2) d found by F13, explaining why we found 
the transits of WASP-57 to occur earlier than predicted. There are 


Table 4. Times of minimum light and their residuals versus the ephemeris 
derived in this work. All but one of the timings were derived in the current 
work, from the source data given in the final column. 


Time of min. 

Error 

Cycle 

Residual 

Source 


(BJD/TDB) 

(d) 

no. 

(d) 




2454602.18313 

0.00317 

-513.0 

-0.00075 

This 

work 

(WASP 2008) 

2454900.26529 

0.00306 

-408.0 

-0.00504 

This 

work 

(WASP 2009) 

2455351.65767 

0.00255 

-249.0 

-0.00071 

This 

work 

(WASP 2010) 

2455686.65086 

0.00055 

-131.0 

0.00009 

This 

work 

(TRAPPIST) 

2455723.55722 

0.00025 

-118.0 

0.00051 

This 

work 

(Euler) 

2455723.55606 

0.00028 

-118.0 

-0.00065 

This 

work 

(TRAPPIST) 

2456001.76950 

0.00035 

-20.0 

-0.00123 

This 

work 

(TRAPPIST) 

2456018.80454 

0.00017 

-14.0 

0.00030 

This 

work 

(Euler) 

2456018.80315 

0.00042 

-14.0 

-0.00109 

This 

work 

(TRAPPIST) 

2456058.54852 

0.00143 

0.0 

-0.00058 

This 

work 

(BUSCA u) 

2456058.54891 

0.00039 

0.0 

-0.00019 

This 

work 

(BUSCA g) 

2456058.54959 

0.00023 

0.0 

0.00049 

This 

work 

(BUSCA r) 

2456058.54913 

0.00042 

0.0 

0.00003 

This 

work 

(BUSCA 2) 

2456728.54099 

0.00171 

236.0 

0.00711 

Dittler (TRESCA) 

2456796.66754 

0.00019 

260.0 

-0.00038 

This 

work 

(DEOSC) 

2456802.35000 

0.00710 

262.0 

0.00424 

This 

work 

(Lopresti) 

2456819.38110 

0.00100 

268.0 

0.00183 

This 

work 

(Corfini) 

2456819.36760 

0.00510 

268.0 

-0.01167 

This 

work 

(Marchini) 

2456833.57422 

0.00017 

273.0 

0.00036 

This 

work 

(DEOSC) 

2456833.57442 

0.00034 

273.0 

0.00056 

This 

work 

(GROND g) 

2456833.57296 

0.00021 

273.0 

-0.00090 

This 

work 

(GROND r) 

2456833.57389 

0.00022 

273.0 

0.00003 

This 

work 

(GROND i) 

2456833.57460 

0.00039 

273.0 

0.00074 

This 

work 

(GROND 2) 

2456836.41480 

0.00290 

274.0 

0.00202 

This 

work 

(LBarbieri) 

2456836.41880 

0.00250 

274.0 

0.00602 

This 

work 

(Corfini) 

2456836.41490 

0.00190 

274.0 

0.00212 

This 

work 

(Lopresti) 

2456836.41960 

0.00310 

274.0 

0.00682 

This 

work 

(Lopresti) 

2456836.41840 

0.00290 

274.0 

0.00562 

This 

work 

(Marino) 

2457171.40240 

0.00350 

392.0 

-0.00277 

This 

work 

(Lopresti) 


several plausible reasons for such a discrepancy to occur, but we 
are not in a position to choose between them. 

The X? of the fit is 1.99, and we interpret this as an indication 
that the uncertainty estimates for the timings are too small. Fig.|4] 
shows the residuals of the times of mid-transit versus the ephemeris 
given above. There is no sign of long-term transit timing variations. 
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Table 5. Parameters of the fit to the light curves of WASP-57 from the JKTEBOP analysis (top). The final parameters are given in bold and the parameters 
found by F13 are given below this. Quantities without quoted uncertainties were not given by F13 but have been calculated from other parameters which were. 


Source 

PA ± Pb 

k 

in 

PA 

Pb 

Euler (2011/06/10) 

0.1198 ±0.0084 

0 . 1 139 ±0.0027 

87.04 ± 1.10 

0.1075 ± 0.0073 

0.01225 ±0.00109 

TRAPPIST (2011/06/10) 

0.1065 ±0.0065 

0.1087 ±0.0027 

89.94 ±1.36 

0.0961 ±0.0055 

0.01044 ±0.00073 

Euler (2012/04/01) 

0.1256 ±0.0066 

0.1190 ±0.0029 

86.56 ±0.73 

0.1123 ±0.0057 

0.01335 ± 0.00099 

TRAPPIST (2012/03/15) 

0.1436 ±0.0095 

0.1263 ±0.0026 

85.28 ±0.84 

0.1275 ±0.0082 

0.01610 ±0.00136 

TRAPPIST (2012/04/01) 

0.1405 ±0.0153 

0.1179 ±0.0042 

85.16 ± 1.20 

0.1257 ±0.0140 

0.01482 ±0.00194 

BUSCA g (2012/05/10) 

0.1370±0.0131 

0.1188 ±0.0056 

85.53 ±1.12 

0.1224 ±0.0113 

0.01455 ±0.00170 

BUSCA r (2012/05/10) 

0.1340 ±0.0057 

0.1186±0.0015 

85.70 ±0.52 

0.1198 ± 0.0050 

0.01421 ±0.00074 

BUSCA ^ (2012/05/10) 

0.1328 ±0.0098 

0.1230±0.0018 

85.69 ±0.84 

0.1182 ±0.0086 

0.01454 ±0.00116 

DEOSC (2014/05/18) 

0.1309 ± 0.0057 

0.1173 ±0.0018 

86.04 ±0.52 

0.1171 ±0.0049 

0.01374 ± 0.00077 

DEOSC (2014/06/24) 

0.1306 ±0.0066 

0.1166 ±0.0017 

86.05 ±0.63 

0.1170 ±0.0057 

0.01364 ±0.00086 

Final results 

0.1278 ±0.0033 

0.1182 ±0.0013 

86.05 ± 0.28 

0.1143 ±0.0029 

0.01331 ± 0.00051 

P13 

0.1122 

0.1127 ±0.0006 

88.0t°;i 

0.1008 

0.01135 


4 LIGHT CURVE ANALYSIS 

Eight of our light curves cover a full transit at high photometric pre¬ 
cision. The Euler and one of the two TRAPPIST light curves from 
P13 also satisfy this criterion. Each of thes e ten datasets was mod- 
elled separately using the JKTEBOlQcode dSouthworth et al. |2004l) 
and the Homogeneous Studies methodology ( ISouthworthl2012L and 
references therein). We did not subject those light curves with only 
partial coverage of a transit to this analysis, because the param¬ 
eters derived from partial light curves are highly uncertain - so 
have little effect on the final results - and are often unreliable (e.g. 
iGibson et alj200^ . 

The JKTEBOP model is based on the fractional radii of the 
star and the planet (ta and rp), which are the ratios between the 
true radii and the semimajor axis (ta.u = The parameters 

of the fit to each light curve were the sum and ratio of the frac¬ 
tional radii (ta + fb and fc = -^j, the orbital inclination (i), limb 
darkening coefficients, and the time of mid-transit. We assumed an 
orbital eccentricity of z ero, based on the finding by F13 that the 
ILucv & SweenevI dl97lh test yielded a 100% probability that the 
orbit was circular. We fixed the orbital period to the value found in 
Section[3 We also fitted for the coefficients of a first-order poly¬ 
nomi al relating differential magnitude and time llSouthworth et al.l 
|2014) . in order to allow for any errors in flux normalisation which 
change with time or airmass. The TRAPPIST light curve obtained 
on the night of 2012/04/01 was split into two sequences by a merid¬ 
ian flip. This was accounted for by modelling both sequences to¬ 
gether but specifying a separate polynomial (of order 1) for each 
sequence. 

Limb darkening (L D) was incorporated using each of five laws 
(see ISouthworthll2008r) . with the linear coefficients either fixed at 
theoretically predicted valued or included as fitted parameters. 
We did not calculate fits for both LD coefficients in t he four two- 
coeffi cient laws as they are very strongly correlated dCarter et al.l 
I2OO8I) . The nonlinear coefficients were instead perturbed by ±0.1 
on a flat distribution during the error analysis simulations, to ac¬ 
count for the uncertainty in theoretical LD coefficients. 

® JKTEBOP is written in FORTRAN77 and the source code is available at 
http://www.astro.keele.ac.uk/jkt/codes/jktebop.html 
® Theoretical LD coefficients were obtained by bilinear interpola¬ 
tion in Tefl and logg using the JKTLD code available from: 
http://www.astro.keele.ac.uk/jkt/codes/jktId.html 


Error estimates for the fitted parameters were obtained in four 
steps. Steps 1 and 2 were resi dual-permutation and Monte Carlo 
simulations dSouthworthlbOOSl) . and the larger of the two alterna¬ 
tives was retained for each fitted parameter. Eor step 3 we ran so¬ 
lutions using the five different LD laws, and increased the errorbar 
for each parameter to account for any disagreement between these 
five solutions. For step 4 we calculated the weighted mean of each 
photometric parameter using the values found separately from each 
light curve. This final step is a powerful external check on the re¬ 
liability and mutual agreement between different datasets, as any 
discrepancies are obvious and quantifiable. 

For all ten light curves we found that it was possible to fit for 
one of the two LD coefficients: reasonable values for the coeffi¬ 
cients were obtained as well as a slightly smaller xS compared to 
fits with both LD coefficients fixed. We therefore adopt these re¬ 
sults, which are summarised in Table|5] Detailed tables of results 
for each light curve are available in the online-only Appendix. The 
best fits are plotted in Fig.|5] 

We find that the results from the different light curves are not 
in perfect agreement, with a xS of 2.0 for ta ± fb and ta, 2.8 for 
Pb and 3.7 for k versus the weighted mean of their values. This is 
due primarily to the TRAPPIST light curve from FI3, which has a 
very small pa and high i compared to the other datasets (TablefSt. A 
degen erac y between these pa rameters is common (e.g. lCarter et al.l 
I2OO8I and Isouthworthi l2008h and arises because these two values 
together specify the observed transit duration, a quantity which is 
well-determined by high-quality light curves. If we adopt instead 
the results from fitting this light curve with both LD coefficients 
fixed, the agreement becomes much better: k has xi ~ 2.5, the 
other four parameters in Table[5]all have Xv < 0.9, and all photo¬ 
metric parameters change by less than their Icr errorbars. We have, 
however, chosen not to take this step for two reasons. Firstly, theo¬ 
retical LD coefficients are not perfect - if they were then different 
sources would give exactly the same values - and none are avail¬ 
able calculated specifically for the TRAPPIST I+z filter. Secondly, 
taking an alternative approach for a discrepant dataset raises the 
possibility of causing an underestimate of the true uncertainties in 
the measured quantities. We have therefore retained the discrepant 
values when calculating the weighted means of the photometric pa¬ 
rameters, and have inflated the errorbars on the weighted means by 
in order to account for the discrepancy. 

The photometric parameters found by FI3 differ significantly 
from our results, by A.la for pa and 3.8 (j for Pb (using our errorbars 
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Figure 5. Phased light curves of WASP-57 compared to the JKTEBOP best fits (left) and the residuals of the fits (right). Labels give the source and passband 
for each dataset. The polynomial baseline functions have been removed from the data before plotting. Only light curves with full coverage of a transit were 
included in this analysis. 


to calculate the a values as error estimates were not provided by 
F13 for these two quantities). This is due to the dependence of the 
F13 solution on only three transit light curves (two complete and 
one only partially covering a transit), all modelled simultaneously, 
of which one was the TRAPPIST dataset we find to be discrepant. 
The value of k found by F13 (0.1127 ± 0.0006) is 3.8cr smaller 
than ours, and is evidence that the error estimates quoted by F13 
are too small (see lSouthworthl2012l and references therein for other 
examples). An alternative explanation is the presence of starspots, 
which is plausible for a star of this temperature. However, no traces 
of spot occultations are seen in our light curves and no rotational 
modulation is seen in the long-term WASP light curves. 


5 PHYSICAL PROPERTIES 


We measured the physical properties of the WASP-57 system us¬ 
ing the results from Section|4] five grids of pred i ctions from the- 


oretical models of stellar evolution llciaretl 20041; Demaraue et al. 

2004 

;|Pietrinfemi et alJl2004l:IVandenBere et alj|2006l: Dotter et al. 

2008 

), and the host-star spectroscopic properties. Theoretical mod- 


els provide an additional constraint on the stellar properties, needed 
because the system properties cannot be obtained from only mea¬ 
sured quantities. The spectroscopic properties were obtained by 
F13 and comprise effective temperature (Teu = 5600 ± 100 K), 
metallicity ([Fe/H] = —0.25 ± 0.10) and velocity amplitude 
{ Ka ~ 100 ± 7m s~ ^). We used the physical constants tabulated 
bv ISouthworthI (1201 ih . 

We first estimated the velocity amplitude of the planet, Kh, 
which was used along with the measured ta, fu, i and Ka to de- 
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Table 6. Derived physical properties of WASP-57. The values found by F13 are given for comparison. 


Quantity 

Symbol 

Unit 

This work 

F13 

Stellar mass 

Ma 

M0 

0.886 ± 0.061 ±0.028 

0.954 ±0.028 

Stellar radius 

Ra 

R© 

0.927 ± 0.031 ±0.010 

0.836+°;^ 

Stellar surface gravity 

loggA 

cgs 

4.452 ± 0.024 ±0.005 

4.574t“;°°s 

Stellar density 

PA 

p© 

1.113 ±0.085 

±.uoo_Q Qg3 

Planet mass 

Mb 

Mjup 

0.644 ± 0.060 ±0.014 

U.O/Z_g Q4g 

Planet radius 

Rb 

R-Jup 

1.050 ± 0.052 ±0.011 


Planet surface gravity 

9b 

ms“2 

14.5 ± 1.5 

18.31=;® 

Planet density 

Ph 

Pjup 

0.521 ± 0.072 ±0.006 

0.8731®;°™ 

Equilibrium temperature 

T' 

eq 

K 

1338 ± 29 

125ll=i 

Safronov number 

0 


0.0522 ± 0.0045 ±0.0006 


Orbital semimajor axis 

a 

au 

0.03769 ± 0.00088 ± 0.00040 

0.0386 ±0.0004 



o PMI'WI llVl“n l ; i| W1/V WmW " Vr ’1/ 

2 3 5 ,Q 20 30 50 


Period (d) 

Figure 6. Scargle periodograms of the SuperWASP data from the 2008, 2009 and 2010 seasons. The orbital period of the system is shown with a downward¬ 
pointing arrow. The stellar rotational period inferred by F13 from its projected rotational velocity is shown by the horizontal errorbar. 


termine the physical properties of the system (ISouthworthll2009h . 
The estimate of Kb was then iterated to find the best match be¬ 
tween the measured ta and the calculated —, and the observed 
Tefi and that predicted by a theoretical model for the obtained stel¬ 
lar mass, radius and [Fe/H]. This was done for a grid of ages from 
the zero-age main sequence to beyond the terminal-age main se¬ 
quence for the star, in 0.01 Gyr increments, and the overall best 
Kb was adopted. The statistical errors in the input quantities were 
propagated to the output quantities by a perturbation approach. 

We ran the above analysis for each of the five sets of theo¬ 
retical model predictions, yielding five different estimates of each 
output quantity. These were transformed into a single final result for 
each parameter by taking the unweighted mean of the five estimates 
and their statistical errors, plus an accompanying systematic error 
which gives the largest difference between the mean and individual 
values. The final results of this process are a set of physical prop¬ 
erties for the WASP-57 system, each with a statistical error and a 
systematic error. The stellar density, planetary surface gravity and 
planetary equilibrium temperatur es can be calculated without re¬ 
sorting to theoretical predictions dS eager & Mallen-Omelas|[200^ : 
ISouthworth et alfcOOTl i lSouthworthllTOlOh . so do not have an asso¬ 

ciated systematic error. 

ISodenheimer et aO ( l2003l . their tables 1 and 2) provided pre¬ 
dicted radii for planets of mass 0.69 Mjup and = 1000 K and 
1500 K. Their radii are 1.01-1.10 Rjup with, and 1.03-1.13 Rjup 


without, a core or additional kinetic heating of the planetary inte¬ 
rior. Both are in very good agreement with the radius of 1.05 ± 
0.05 Rjup we find for WASP-57 b. 

Table|^ contains our measurements of the physical properties 
of the WASP-57 system. Compared to F13, we find a less mas¬ 
sive but larger star. As planetary properties are measured relative 
to those of their parent star, the planet is similarly affected. The 
measured planetary density is 3.5 ct lower, at 0.521 ± 0.072 pjup 
compared to the value of O.STSlQg^f pjup found by F13. Our re¬ 
sults are based on a much more extensive set of photometric data 
so are to be preferred to previous measurements, even though the 
errorbars have not changed by much. A significant advance in our 
understanding of the WASP-57 system could be achieved by ob¬ 
taining further spectroscopy of the host star, from which more pre¬ 
cise values for Teu, [Fe/H] and Ka could be measured. This is of 
particular interest because of its metal-poor nature, whose effect on 
the i ncidence of different types of planets is curren tly under discus¬ 
sion teuchhave et alJlToI^Jwang & Fischej2015h . 

The age of the system is unconstrained in our analysis above, 
as often occurs when the host star is significantly less massive 
than 1 Mq. F13 inferred age estimates of ^2 Gyr for WASP-57 A 
from its photospheric lithium abundance, and ~1.9jl^'2 frorn 
gyrochronological arguments and its rotation period derived from 
its radius and projected rotational velocity. The Teu of the star is 
within the regime where starspots are common so we have checked 
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Table 7. Values of and for each of the light curves. The errorhars exclude all common sources of uncertainty in rj-, and R\^ so should only be used to 
compare different values of r’i,(A). The final column gives the size of the errorbar on R\^ in atmospheric scale heights. 


Instrument 

Passband 

'^cen 

(nm) 

FWHM 

(nm) 

rh 

Rb 

(Rjup) 

a{H) 

Euler 

Gunn r 

660.0 

100.0 

0.01362 ±0.00004 

0.867 ±0.038 

8.2 

BUSCA 

Gunn u 

350.0 

68.0 

0.01099 ±0.00048 

1.061 ±0.008 

1.8 

BUSCA 

Gunn g 

495.5 

99.5 

0.01345 ±0.00011 

1.063 ±0.005 

1.2 

BUSCA 

Gunn r 

663.0 

105.0 

0.01348 ±0.00007 

1.100 ±0.011 

2.3 

BUSCA 

Gunn 2 

910.0 

90.0 

0.01395 ±0.00013 

1.057 ±0.017 

3.6 

DFOSC 

Bessell R 

648.9 

164.7 

0.01343 ±0.00003 

1.060 ±0.009 

2.0 

GROND 

Gunn g 

477.0 

137.9 

0.01340 ±0.00021 

1.033 ±0.010 

2.2 

GROND 

Gunn r 

623.1 

138.2 

0.01344 ±0.00011 

1.036 ±0.017 

3.1 

GROND 

Gunn i 

762.5 

153.5 

0.01310 ±0.00013 

1.059 ±0.003 

0.6 

GROND 

Gunn z 

913.4 

137.0 

0.01314 ±0.00022 

1.074 ±0.004 

0.8 


if it is possible to precisely determine its rotation period from spot- 
induced modulation. A Lomb-Scargle periodogram was calculated 
for each of the three seasons of SuperWASP data and can be seen 
in Fig.|^ There are no strong peaks in the period interval of inter¬ 
est (5-30 d), and no moderately strong peaks present at the same 
period in all three seasons. We conclude that the rotational mod¬ 
ulation of the star is below the level of detection with the current 
data. 


5.1 Comparison with theoretical models of giant planets 

F13 found that the measured mass and density of WASP-57 b im¬ 
plied the presence of a heavy-element core of mass ro ughly 50 Mf^, 
via a comparison to the theoretical predictions of IFortnev et al.l 
ll2007l) . This rather large core mass is surprising given the sig¬ 
nificantly sub-solar metal abundance of the host star ([Fe/Fl] = 
—0.25 ± 0.10). As we have found a significantly lower density for 
the planet (smaller by 40% or 3.5 (j) it is germane to reconsider this 
conclusion. We have therefore compared our new mass and radius 
measurements with predictions based on three batches of theoreti¬ 
cal models. _ 

ISaraffe et alj (120081 their table 4) find planetary radii of 0.97- 
1.06 Rjup for planets of mass 0.5-1.0 Mjup and age 0.5-5 Gyr, 
again in accord with our results. These values are for a heavy- 
element fraction of Z = 0.02, and larger fractions result in pro¬ 
gressively smaller radii and thus poorer agreement with our radius 
measurement. 

Finally, the pro perties of WASP-57 b match the predictions of 
iFortnev et al.l l l2007l their fig. 6) for a 25 M® heavy-element core. 
The difference in radius between models with and without this 
core are only 0.05 Rjup for a 1 Mjup planet and 0.18 Rjup for a 
0.3 Mjup planet, so are of a comparable size to the uncertainty in 
the radius of WASP-57 b. Our measured properties for this planet 
therefore do not provide significant support for a high metallicity 
or the presence of a heavy-element core. 


6 VARIATION OF RADIUS WITH WAVELENGTH 

Two of our datasets include observations in four passbands si¬ 
multaneously (ugrz for BUSCA and griz for GROND), whereas 
we have r- or i?-band photometry from four different sources 
(DFOSC, EulerCam, BUSCA and GROND), so it is relevant to 
search for possible changes in the measured radius of the planet as 


a function of wavelength. Such anal yses are the photometric e quiv- 
alent of transmission spectroscopy dSeager & Sasselo^|2000|) and 
have been pioneered at op tical wavelengths by Sing et al. 1201 1 ab ¬ 
ide Mooii et alj (12012h and I Southworth et"^ ( 1201 2 I 1 . 

Changes in the radius measured from planetary transits, as 
a function of wavelength, are predicted to occur due to opacity 
variations which affect the height at which the atmosphere trans¬ 
mits light coming from the parent star in the direction of the ob¬ 
server. At blue wavelengths a greater atmospheric opacity due to 
Rayleigh and Mie scattering leads to a higher maximum depth at 
which starlight is trans mitted, causing an increase in the mea sured 
radius of the planet (e.g. lPont et al.l200^ : lNikolov et alj20T^ . En¬ 
hanced opacity also le ads to signatures of sodium and potassium 
at optical wavelengths dEortnev et alj|2008h. although only narr ow 
absorption cores have been detected so far dNikolov et al.ll201^ . 

The fundamental observable in this work is the transit depth, 
represented in our notation by the ratio of the radii k or the frac¬ 
tional planetary radius ru. The parameter directly comparable to 
theoretical predictions is the true planetary radius, Rh- The param- 
eter ru is correlat ed with other photometric parameters (see e.g. 
ISouthwortiilbOOSh . and its transformation into Rh requires other 
parameters which have uncertainty but are common to all photo¬ 
metric passbands (the error budgets calculated in the previous sec¬ 
tion show that these are ta, Tefi and [Fe/H]). 

We removed these two effects in order to determine Rh values 
with relative errorbars. We did this by refitting the light curves with 
all parameters fixed except k. To, the linear LD coefficient for the 
quadratic LD law, and the coefficients of the polynomials between 
differential magnitude and time. We then transformed the result¬ 
ing rb values into Rh using a fixed orbital semimajor axis, a. This 
yielded a set of Rh values and errorbars which are directly com¬ 
parable to each other. The uncertainties in rb were measured using 
1000 Monte Carlo simulations each. 

The data included in this analysis were the two Euler light 
curves (modelled simultaneously), the two DFOSC light curves 
(modelled simultaneously), the BUSCA ugrz and the GROND 
griz data. We did not include the TRAPPIST light curves because 
the very wide passbands (I+z or blue-blocking filters) yield mini¬ 
mal spectral resolution. Although the GROND data only partially 
cover a transit, which precluded their use in Section|4] they give 
reliable results here because ta and i were fixed during the fitting 
process. 

Fig.|3 shows the resulting values of Rb as a function of the 
central wavelength of the passbands used. The FWHMs of the pass- 
bands are shown for reference using horizontal lines. The originat- 
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Figure 7. Measured planetary radius (iJb) as a function of the central wave¬ 
length of the passbands used for the different light curves. The datapoints 
show the iJb measured from each light curve. The vertical errorbars show 
the relative uncertainty in fJi, (i.e. neglecting the common sources of error) 
and the horizontal en'orbars indicate the FWHM of the passband. The dat¬ 
apoints are colour-coded according to passband, and the passbands are la¬ 
belled at the top of the figure. The symbol types are filled circles (BUSCA), 
open circles (GROND), upward-pointing arrow (DFOSC) and downward¬ 
pointing arrow (Euler). On the right of the plot we show the value of iJr, 
measured in Section|4]and the size of ten atmospheric pressure scale heights 
(lOiT). The grey lines through the empirical datapoints show theoretical 
predictions for a transmission spectrum of a gas-giant planet of solar chem¬ 
ical composition from Madhusudhan (priv. comm.). The darker-grey line 
includes features due to Na and K whereas the lighter grey line also in¬ 
cludes TiO opacity. 

ing rb values and passband characteristics are collected in Tabled 
Two conclusions are immediately apparent from this figure. Firstly, 
the planetary radius in the w-band is very small and very uncertain. 
Secondly, the grRiz results are consistent with no variation of i?B 
with wavelength. 

The M-band light curve shows a small transit depth and a high 
scatter (see Fig.|7}, which causes the anomalous rb measurement 
in this band. The discrepancy relative to the overall rb value ob¬ 
tained in Section|4]is highly significant at 4.8cr, and corresponds to 
approximately SfiTT. H is the pressure scale height, and in the case 
of WASP-57 b is 334 ± 35 km (0.0047± 0.0005 Rjup). A variation 
in i?b of the size of 39iT is difficult to explain, and is not believ¬ 
able unless confirmed by additional data of much higher quality. 
The two 2 -band light curves also show a clear disagreement of size 
3.2cr. 

Fig.[7]also shows two theoretical transmission spectra calcu¬ 
lated under different assumption s by Madhusudhan (priv. comm .) 
using the atmosphere code of iMadhusudhan & SeageJ j2009t) . 
These predictions are for a gas-giant planet of radius 1.25 Rjup and 
surface gravity 25 ms“^ so have been scaled to match the smaller 
radius and lower gravity of WASP-57 b. Our finding of a small u- 
band rb is not consistent with these theoretical predictions. 

A close-up of the main part of Fig.[7] is shown in Fig.[8] 
Passband-averaged values for the transmission spectra are shown 



Figure 8. Close-up of the main part of Fig.|3showing the 7?b measurements 
and theoretical transmission spectra. The u-band result is off the plot. Its 
central wavelength is indicated with a downward-pointing arrow. The size 
of one atmopsheric scale height is indicated to the right of the plot. The 
black circles are the values of passband averages of the two transmission 
spectra, and are shown at the central wavelengths of the relevant passbands. 

with black filled circles, and differ by up to 2H. Th e Rayleigh 
scatteri ng slope could be significantly greater than this: ISing et al.l 
l l2011bl) found that Rayleigh scattering caused the measured radius 
of HD 189733 b to be larger by 577 at 400 nm than at 900 nm. 

The relative uncertainties in our measured radii for WASP- 
57 b are below 177 for two, and below 277 for five, of the ten light 
curves (see Table|7j. We are therefore sensitive to radius variations 
at the level of 177, which is smaller than both the difference be¬ 
tween the two theoretical transmission spectra and the size of the 
Rayleigh scattering slope detected for HD 189733 b. Our data are 
therefore sensitive, in principle, to the atmospheric properties of 
WASP-57 b. However, in practise, our measurements are insuffi¬ 
cient for studying the atmosphere of this planet due to the anoma¬ 
lous result for the «-band and the scatter of the radius measure¬ 
ments in the i- and 2 -bands. The situation could be improved by ob¬ 
taining data in narrower passbands (i.e. higher spectral resolution), 
and with repeated observations over the full optical wavelength 
range. Particular attention should be paid to the n-band, which is an 
important discriminant between the two transmission spectra and 
also enhances sensitivity to the Rayleigh scattering slope. 


7 SUMMARY AND CONCLUSIONS 

WASP-57 b is a relatively low-mass hot Jupiter orbiting a cool star. 
Amateur astronomers first noticed that its transits were occurring 
earlier than predicted, a finding subsequently confirmed by obser¬ 
vations from professional facilities. We have presented ten transit 
light curves from amateur astronomers, plus 13 obtained using pro¬ 
fessional telescopes of which seven predate the discovery of inac¬ 
curacy in the orbital ephemeris of the system. We have determined 
a revised orbital ephemeris which differs by 24a from the orbital 
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period in the discovery paper, and can be used to predict transits 
to a precision of less than 1 minute until the year 2170. We also 
obtained high-resolution Lucky Imaging observations, which show 
no evidence for nearby companions whose flux might have contam¬ 
inated our light curves. 

We have used these and previously published data to redeter¬ 
mine the physical properties of the WASP-57 system, finding that 
both the planet and its host star are larger and less massive than 
previously thought. A comparison of our new results for WASP- 
57 b to theoretical predictions for the properties of gaseous planets 
reveals a good agreement with models lacking a core or additional 
heat sources. This disagrees with the core mass of 50 M® postu¬ 
lated by FI3, but is in accord with expectations for a planet which 
formed around a star of significantly sub-solar metal abundance. 

We observed two of the transits of WASP-57 using two 2.2 m 
telescopes equipped with simultaneous multi-band imaging instru¬ 
ments: GROND {griz passbands) and BUSCA {ugrz passbands). 
These data are well-suited to investigating the possible wavelength- 
dependence of the planet’s measured radius due to effects such as 
Rayleigh and Mie scattering, and atomic and molecular absorption. 
Whilst the radii in the g and r!R bands are in generally good agree¬ 
ment, the i-band measurement is slightly smaller than expected and 
the two z-band measurements are discrepant by 3 (t. The u-band 
radius is crucial for measuring the Rayleigh scatte ring slope, as 
well a s separating the pM and pL classes proposed bv lportnev et aU 
( l2008h . Our measurement is 50" below theoretical predictions, and 
the size of the discrepancy is inexplicable using current theoretical 
transmission spectra. This result is almost certainly spurious, and 
can plausibly be blamed on the strong absorption by Earth’s atmo¬ 
sphere at blue-optical wavelengths plus the faintness of the host star 
in this passband. 


7.1 Future opportunities for transmission photometry 

Successful detections of radius variations in optical transmis¬ 
sion photo metry have recently b e en announced for the TEPs 
GJ 3470 b dNascimbeni et Zll2013l: [Biddle et alj|20l4. Oatar-2 b 


dlVtancini et alj 2014 ) and WASP-103 b l lSouthworth et alJl201^ . 


Only one of these studies presented data obtained shortward of the 
Balmer jump, which is an important but observationally difficult 
wavelength interval (see Fig.[^. 

Transit light curves in the u and U bands have previ¬ 
ously been presented for several TEPs, and have shown plane¬ 
tary radii_either_con^tent_with other optical passbands (WASP- 
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Copperwheat et al.ll2013l : TrES-3. [Turner et al.ll2013l: WASP- 


Bento et al.ll20l4 WA SP-39 and WASP-43. iRicci et alj[20ia: 
XO-2. IZellem et ail 2 OI 5 I 1 or somewhat lar ger than other optical 
passbands (HAT-P-5, ^mhworth et alj2012L J. Dittman 2012, priv. 


comm.; GJ 3470. iNascimbeni et al. 2013h . A universal feature of 
these studies is the reliance on either a single ri-band transit light 
curve, which yields large uncertainties on the measured planetary 
radius, or the use of data not obtained simultaneously in multiple 
passbands, so the results are hostage to temporal changes such as 
induced by magnetic activity in the host stars. Most transmission 
photometry studies also suffer from the use of wide passbands, 
which are insen sitive to spectral fea tures other than broad contin¬ 
uum slopes (see lNikolov et al.|l2013h . 

Whilst suffering from a lower spectral resolution, transmission 
photometry has several advantages over transmission spectroscopy. 
These include being able to observe over a wide wavelength inter¬ 
val without being subject to second-order contamination, the ability 
to use comparison stars more distant from the planet host star, and 


the option to use tele scope defocussing techniques to avoid sys¬ 
tematic noise (but see lBurton et alJl2015L for a counter-example). 
Smaller telescopes can be used, making it easier in particular to ob¬ 
serve multi ple transits and thus demonstrate the repeatability of the 
experiment teean et alj2013l;lGibs^l2014h . 

We therefore advocate studies based on observations of mul¬ 
tiple transits, obtained simultaneously through many intermediate 
or narrow passbands. These passbands should be well-defined by 
interference filters, thus avoiding compromises such as the vari¬ 
able red edge of the z filter due to its reliance on the quan tum 
efficiency curve of the CCD used (e.g. lEukugita et al.l 1 19961) or 
the red leak in some u and U filters which is capable of causing 
spurio us results for optically-blue objects (e.g. lGuhathakurta et al.l 
Il998h . With a sufficient number of passbands, it should be possible 
to achieve low-resolution spectroscopy of the atmospheres of ex¬ 
trasolar planets through the full optical wavelength range using the 
transmission-photometry approach. 
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